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Interorbital pair scattering in clean and impure superconductors 
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The calculations of the local and global properties of two band superconductors have been pre- 
sented with particular attention to the role of the inter-orbital scattering of pairs. The properties 
of such superconductors are very different from a single band or typical two band systems with 
dominant intra-band pairing interactions. The role of Van Hove singularity in one of the bands on 
the properties of intra-band clean superconductor has been discussed. It leads to marked increase of 
superconducting transition temperature in the weak coupling limit. We study the inhomogeneous 
systems in which the characteristics change from place to place by solving the Bogolubov-de Gennes 
equations for small clusters. The suppression of the superconducting order parameter by the single 
impurity scattering the fermions between bands is contrasted with that due to intra-band impurity 
scattering. The results obtained for impure systems have been shown as a maps of local density of 
states, the order parameter and gap function. They can be directly compared with STM spectra of 
the real material. 

PACS numbers: 71.10.Fd; 74.20.-z; 74.81.-g 



I. INTRODUCTION 

Already in the fifties and sixties the main properties of 
the two band superconductors have been clarifiedii^i^ii^. 
At that time, however, the existing materials did not 
show clear evidence of two band behavior. The experi- 
mental situation has changed with the discovery of the 
high temperature superconducting oxides^ and even more 
with subsequent discoveries of strontium ruthenate^, 
magnesium diboridc^ and iron pnictides^i^. Even though 
all of these systems have a number of bands in the vicin- 
ity of the Fermi energy their presence shows up in quite 
a different way. 

Magnesium diboride clearly shows two different gaps 
of the same symmetrj*iSiii. In strontium ruthenate 
the three band model seems to be necessary to ex- 
plain its puzzling propertie o^^'^'^ . The model of super- 
conductivity in the iron pnictides is a matter of ongo- 
ing debatei^ii^ii^iii. The iron pnictides possess a large 
number of bands around the Fermi energy and few of 
them seem to play important role in the superconduct- 
ing statoi^iiSii^. Two band model has been proposed as 
minimal model of these superconductor^^ . 

With two bands near the Fermi energy one generally 
expects formation of intra-band and inter-band pairs. In 
the later case the pairs have in general non-zero center of 
mass momentunii^^. The simpler case^ of superconductiv- 
ity with intraband pairs, which can be scattered between 
two bands seem to be relevant in modelling of pnictides. 
Indeed, there are strong theoretical^^i^i arguments that 
inter-band interactions may be important in these sys- 
tem. These findings make pnictides different from MgB2, 
in which main coupling mechanism is intraband^^. Thus 
the detailed study of inter-orbital^- pair scattering mech- 
anism of superconductivity is timely and of importance. 
The issue has recently been discussed in connection with 
both cuprate^? and pnictide superconductors^^. 



In this paper we are mainly interested in the proper- 
ties of inter-orbital only mechanism of superconductiv- 
ity. We shall study both clean homogeneous and impure 
models. Without loss of generality we shall denote two 
orbitals as 1 and 2. The interaction Un {U22) is responsi- 
ble for superconducting instability inside a band formed 
by orbitals 1 (2), while U12 promotes the scattering of 
superconducting pair between orbitals 1 and 2. The im- 
purity scattering potential is assumed in general form 

^mp'^iAcr'^iAV- It scatters electrons from site i, orbital 
A' into orbital A of the same site. If A = A' we call such 
impurities intra-band, if A 7^ A' inter-band. 

The organisation of the rest of the paper is as follows. 
Section 2 presents the general Hamiltonian of the two 
orbital model and the Bogolubov - de Gennes (BdG) ap- 
proach used to solve it. The homogeneous superconduc- 
tors are discussed in Section 3, where we study inter alia 
the effect of Van Hove singularity in the density of states 
in one of the bands on the properties of the superconduc- 
tors with inter-band pair scattering only. The changes 
induced in the superconductor by single intra-band or 
inter-band impurity are discussed in Section 4, while the 
finite concentration of impurities is considered in Section 
5. We end up with the discussion of our results and their 
relevance to most prominent two band superconductors: 
MgB2 and iron pnictides. 



II. HAMILTONIAN FOR THE TWO ORBITAL 
SUPERCONDUCTOR 



We start with general Hamiltonian in a real space de- 
scribing the system with two orbitals. We assume the 
spin-independent effective pairing interaction between 
fermions in various orbital states. The randomness in 
the system is easily incorporated via site dependence of 
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parameters. The Hamiltonian reads 

ij,XX' .(T 
i,X.{T 

i,Ai A2 jAs A4 

where c^g.; ^iAo- are creation and annihilation opera- 
tors of electrons with spin a =1,1 at the lattice site 
fl = i m the orbital A. ex is the electron energy and 
/i is the chemical potential, t^j^ are the hopping inte- 
grals between the same or different orbitals (if A 7^ A'). 
^AiA2A3A4 (?^) denotes interactions, which are attractive 
for t/AiA2A3A4 (j^) < 0. The dependence of the interaction 
parameters on the position fl allows to treat systems with 
inhomogeneous pairing. 

We use standard mean-field decoupling valid for a spin 
singlet superconductor and get the following effective 
Hamiltonian 

ij,XX' ,cr 
i,X,(T 

+ ^ (A,v(rl)c+^c+,^ + /j.c.), (2) 

i,AA' 

where the order parameters Aaa'I*^) are related to the 
pairing correlation functions /aa'(^) CiAj,CiA't > 

through 

^AiA2(''l) = - X! C^AiA2A3A4(''l)/A3A4(n). (3) 
A3 A4 

The local Hartree terms Vxlfl) depend on the number of 
particles at given site nxaifl) =< cl^^axa > 

Vxain) = ;^C/A'AAA'('^l)nA'-.('^l)- (4) 
A' 

We consider here only diagonal correlations 

< C^o-CiA'cr' SxX'Sacr'nXairi). 

The Hamiltonian ([2]) is diagonalised with help of the 
Bogolubov - Valatin transformation's*^ 

c^AT = ("A,.(rl)7i.T - v*xArl)j+^'j , (5) 
^'Ai = Y ("A..(rl)7L/i + t'L(n)7j|) (6) 

leading to the Bogolubov- de Gennes (BdG) equations 
for amplitudes uxvifi), vxv(fi) and eigenenergies E^, 

Ej,A' Kl^^x'Ar,) + Ea' AAA'(^l)«A'.(rl) = 

E^ux^ifi), (7) 
-E,,A'^,f «Av(r,)+EA' A*,,(rl)uv.(rl) = 

E.vxAn), (8) 



where the operator K^j^ reads 

K^^' = {ex-fi + VxAnM.Sxx' + V^X(^-l)6,, - t^/'.(9) 

The pairing parameters Axx'i'H,) and Hartree poten- 
tials Vxifl) are in turn expressed in terms of eigenfunc- 
tions and eigenergies uai/('^), vxi^ifl), E^, as^^ 

Mrl) = E (l"A.(rl)P/. + l«A.(rl)P(l - /.)) , (10) 

fw'in) = ^[uA^(rl)v^,^(rl)(l - /^) 

- ux'u{ri)v*x^{fi)fu]- (11) 

In the above formulae = ^gE^/ksT ^ ^ denotes 
the Fermi-Dirac distribution function of quasi-particles. 
The total number of particle in the given band (to be 
denoted by the same index as the orbital) is given by 
Nx = E.nxirl). 

The local density of states (LDOS) N{n,E) is directly 
accessible in scanning tunneling microscope (STM) mea- 
surements and is proportional to the local conductance 
dl{fi, V) /dV . In the two band system it is a sum of local 
densities of states of the individual bands A^(A, r^, E) 

N{X,n,E) = ^[|u,A(rl)P<5(^-^,) 

+ Kx{fi)\''5{E + E,)]. (12) 

Obviously we have at each site N{ri,E) = Nil^fl^E) -|- 
N{2,fi,E). For a clean system equations (|7I8|) can 
be Fourier transformed and written (in closely analo- 
gous form) in reciprocal space. For the impure systems 
with broken translational symmetry the Bogolubov - de 
Gennes equations (|7l8p are solve self-consistently in real 
space for a small nxm cluster with periodic boundary 
conditions. For a two orbital model the typical size of 
the cluster is 20x30. In the next section we start with 
the comparison of our real space (for small cluster) cal- 
culations with (numerically) exact results obtained in re- 
ciprocal space {i.e. for bulk system). 

III. HOMOGENEOUS SUPERCONDUCTORS 

In this section we shall discuss some properties of ho- 
mogeneous two band superconductors paying special at- 
tention to the comparison of the accuracy of small cluster 
calculations with bulk system. We also consider the ef- 
fect of Van Hove singularity in one of the bands on the 
properties of inter-band pairing superconductivity and 
the role of various inter-band couplings. 

A. Small clusters vs. bulk systems 

We start with the homogeneous system with two or- 
bitals denoted 1 and 2. The superconductor is de- 
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FIG. 1; Energy bands of the two dimensional square lattice 
with parameters 62 — ei = 2t, ti = t,t2 = 2t, t\2 ~ 0.054. The 
chemical potential is located at jj. — and the resulting total 
number of particles n = 1.62. 



scribed by the following set of parameters. The inter- 
band interaction^ has the form of pair scattering only 
C/12 = V\Y22i while two intra-band interactions are \J\\ — 
and J722 — Vtiti- We consider two-dimensional 
square lattice with non-zero hopping integrals between 
the nearest neighbor sites only t\ — and hybridiza- 
tion ti2 = t]^ . We set the direct hoping between orbitals 
no. 1 as our energy unit t\ = t = \. Since we are ignoring 
the possibility of inter-orbital pairs, we use the simpler 
notation Ai — An and A2 — A22. 

Fig. ([T|) shows the single particle energy bands along 
main symmetry direction in the two dimensional Bril- 
louin zone obtained for the following set of parameters 
62 — 61 — 2t, ti = t,t2 = 2t, ti2 = 0.05i. The chemical 
potential fj, = and the total number of carriers n = 1.62. 

Fig. ^ compares the solutions obtained for the bulk 
system with those for small clusters of various size. We 
consider here the bulk data as exact. The accuracy of de- 
termination of the gap parameter for bulk homogeneous 
system (in our case assumed to be of the order of 10~^t) 
is only limited by the time of calculations. We have found 
that at the band center the results obtained for clusters 
with size greater than 400 sites are acceptable. Relative 
changes of the gaps with respect to the bulk values 



^bulky^bulk . ^QQy^ 



are in the range of SAi < 0.15% in the first band and 
slightly greater SA2 < 1.5% for the second band. Well 
inside the bands the spectrum is quasi-continuous and 
the results agree very well with bulk data but near the 
band edges the spectrum of finite clusters is discrete and 
the differences are larger (c.f. Fig. ([2])). 



Inter-orbital pairing only superconductor 
role of Van Hove singularity 



the 



As a general rule one finds that the inter-orbital scat- 
tering U12 plays a minor role in superconductors with 
dominant intra-band interactions^^,. This interaction, 
however couples two bands and may lead to the increase 
of the superconducting transition temperature^!^. The 
situation changes drastically if the inter-band pairing is 
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FIG. 2: The upper panels show the dependence of the A a 
for both bands on the chemical potential /i for bulk system 
(solid curve) and for clusters of size 15x19 (thin line with 
plusses) and 21x27 (line with crosses). The other parameters 
are: 62 - ei = 2t,ti ^ t,t2 = 2t,ti2 = 0.05t, Ui = -3.54, 
U2 = —34, U\2 = —0.54. At the lower panels temperature 
dependence of Aa are shown for ei — /i = 24. The insets show 
the data near T=0 on the expanded scale. 



the only existing interaction. The properties of the su- 
perconductors with dominant inter-band scattering are 
markedly different from those with dominant intra-band 
interactions. In particular'^'*, the superconducting tran- 
sition takes place for arbitrary sign of C/12. The value of 
the gap in the first band is determined by the interaction 
C/12 and the density of states in the second band and vice 
versa, the gap in the second band is proportional to the 
partial density of states (DOS) at the Fermi level in the 
first band. 

This can easily be seen from the general two band BCS 
equations* 

Ai(l + C/iiFi) = -[/12A2F2, 



A2(l + t/22F2) = -U12A1FU 



where 



Fa 



f^'^- tanh 
dENx{E) 



2kBT 



(13) 



(14) 



and N\ (E) denotes single particle density of states in the 
band A. In this discussion we interested in the limit of 
inter-band pair scattering only. For Un = U22 = the 
above equations reduce to 



Ai 



-C/12A2 



dEN2{E)- 



tanh 



2kBT 



VE^+^2 



4 



A2 = -iJisAi 



^'""^ tanh 
dENi(E) = 





2/cbT 



A? 



(15) 



It is clear from equations (|T5|) that the value of A in 
the second band is determined by the density of states 
in the first one and vice versa. It is also obvious that 
the nonzero solutions can be obtained for both signs and 
arbitrary small value of the coupling Ui2. For positive 
value of it the order parameters in the two bands have 
oposite signs, while for negative U12 they are of the same 
sign. The inter-orbital pairing only model has a number 
of unusual features. It has been foun d^^i'^^ that the ratio 
A2/A1 = y/Ni/N2, where A^2(^i) is the density of states 
in band 2(1) at the Fermi level and the superconductiong 
transition temperature of the system is given by the BCS- 
like expression 



tiuic / — 1 
Tc = 1.136- — exp 



(16) 



with Ae// = Ao = y^UyV^. 

It often happens that the Fermi level in superconduc- 
tors lies close to the Van Hove singularity. In layered 
systems with nesting properties of the (quasi - two di- 
mensional) Fermi surface the density of states near the 
Van Hove singularity changes logarithmically 



N{E) = No \n{2W/\E\)e{\E\ - W) 



(17) 



with 2W being the band width and Q{x) the step func- 
tion. It is known that the existence of such singularity 
modifies22. the BCS expression for the transition temper- 
ature ()16p . In particular, in the one band case and the 
weak coupling limit^^ A <C 1 it leads to increase of the 
superconducting transition by changing the effective in- 



teraction: -T^— - 

Here we assume the density of states in the second 
band to be singular N{2,E) = N2\n{2W/\E\) near the 
Fermi energy, while that of the first band flat N{\,E) = 
A'^i. Near Tc equations psp are linearised, we approx- 
imate tanh(.x) = min{x, 1) and find (we use here h — 
ks = 1) 



Ai 



-A2C/i2iV2[l + ln 



2W 



In- 



2Tr_ 



, 2W , ujc 
-{- In In 

uJc 2Tc 



i(ln^)^ 
2 ^ 2Tc' 



A2 - -Ai!7i2iVi[l + ln^]. 



(18) 



The analysis of the above set of equations in the weak 
coupling limit {U12 — > 0) leads to the approximate 
BCS like expression for the superconducting transitions 



temperature with 



<!// 



)^/^ and to 



the modification of the prefactor, which changes from 
1.136wc to 1.136a;c(2M^/'^c)^^^- K is interesting to note 
that up to the prefactor the Van Hove singularity in one 



of the bands increases Tc of intra-band superconductor 
at the very weak coupling only Aq ^ 1. 

In a similar way one can calculate the effect of Van 
Hove singularity on the ratio of the gaps ^ at zero tem- 
perature. One finds 



Ai(0) = 



-A2(0)C/l2iV2[l 

1, A2(0). 
2' 



In 



A2(0) 
2W 



1/, .,2 



A2(0) = -Ai(0)?7i2A^iln 



2ujc 
A^ 



(19) 



Even in the extreme weak coupling limit Aq when 
Tc, Ai,A2 the gap ratio is not given by the ratio 
of the densities of states and depends on Tc and thus 
Aq. The ratio ^ decreases from the value much larger 

than for small Aq to values below for larger 

Aq. However, it is interesting to note that the correct de- 
scription of the intra-band superconductivity requires the 
strong coupling theory^-, even if Aq < 1 and Van Hove 
singularity plays similar role in Eliashberg equations^. 



C. The role of the band couplings 

Before the presentation of the real space local prop- 
erties of the model with general interactions we spent 
here some time on discussing the homogeneous systems 
and the influence of model parameters on the supercon- 
ducting bulk state. In particular, we are interested in the 
dependence of superconducting state on the couplings be- 
tween bands. The hybridization parameter ti2 provides 
single particle coupling and the inter-band pair scattering 
U12 provides the direct two body inter-band interaction. 
The hybridization changes single particle spectrum and 
this influences the superconductivity. 

Fig. ([SI) (left panel) shows the changes of the order pa- 
rameter in the first band due to increase of hybridization 
<i2- The strong decrease of Ai with ti2 results from the 
changes of the single particle spectrum in the first band. 
The second band is essentially decoupled, as U12 = 0. 
One observes strong decrease of the projected density of 
states around the Fermi level in both bands. This is 
illustrated in the figure Left panel shows the to- 

tal (dashed curve) and projected onto orbital 1 and 2 
densities of states in the system without any inter-band 
coupling {ti2 — 0), while in the right panel for strong 
hybridization t\2 — 5i. 

As mentioned, the inter-band interaction alone leads 
to the superconducting instability independently if it is 
repulsive or attractive. It induces gaps in both bands. 
The results are shown in the right panel of Fig. ([3]) . The 
value of the gap in the second band is larger because the 
density of states near Ep in the first band is larger (c./. 
equations ^E^). 

The simultaneous presence of the inter-band (U12) and 
intra-band (here Un only) interactions results in an in- 
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FIG. 3: Increase of the hybridization results in strong decrease 
of Ai for ti2 > 2 (left panel). For a given set of parameters 
(e2 — ei = 2t, ti = t,t2 = 2t, Uu = —3.5t, n = 1.2) it is mainly 
due to diminishing of the density of states near the Fermi 
energy (cf. Fig. Q). Right panel shows the dependence 
of order parameters in both bands on the modulus of the 
intraband interaction U12 in a model with ti2 = and Uu — 
U22 = 0. 





FIG. 4: Changes in the normal density of states (total - 
the dashed curve and projected onto respective orbitals- solid 
curves) for the system as in the Fig. ((Sj with the hybridiza- 
tion ti2 = (left panel) and 5t (right panel). The Fermi 
energy is at E=0. 



crease of the order parameter in the active band, the 
appearance of A in the nonactive band and character- 
istic modifications of the quasiparticle density of states 
as illustrated in Fig. ([5]). The nonzero density of states 
around chemical potential for U12 = is simply a result of 
the absence of couplings between bands (as also ti2 = 0) 
and the lack of pairing in the second band. This presents 
(slightly artificial) case of coexisting normal electrons and 
coherent Cooper pairs in the system. 

This ends up our analysis of the homogeneous two band 





FIG. 5: The dependence of Ai and A2 on U12 in the super- 
conductor with intraband attraction Uu = —3.5t (left panel). 
The right panel shows the energy dependence of the quasipar- 
ticle density of states for few values of the interband interac- 
tion between bare bands (ti2 = 0). The structure in the 
quasiparticle DOS is due to the gap induced in the non-active 
band. 



superconductors. In the next sections we study various 
inhomogeneities starting with single potential scatterers. 



IV. SINGLE POTENTIAL IMPURITY IN A 
CLEAN SUPERCONDUCTOR 



In this section we study a single short-ranged non- 
magnetic impurity embedded in an otherwise clean sys- 
tem. We solve BdG equations ([7]) and ([5]) on a small 
cluster of size L = 13 x 17 with an impurity placed in its 
center. In the two orbital model the impurity may scatter 

electron from a given orbital to the same orbital (intra- 

1 f 2) 

orbital scattering to be denoted V^^^' ) or to other orbital 
(inter-orbital scattering - V^^p) located at the same site. 
The inter-band scatterers were intensively studied^ in 
connection with MgB2. It has been found that Eliash- 
berg theory leads to much slower rate of Tc suppression 
than predicted on the basis of BCS treatment. Here we 
allow for both, the intra-band and inter-band pairing in- 
teractions and compare the Tc changes induced by two 
types of impurities. Instead of strong coupling Eliash- 
berg approach we are using Bogolubov-dc Gennes ap- 
proach which allows for the distortion of the wave func- 
tion around impurity and is more suitable to treat inho- 
mogeneous superconductors than either BCS or Eliash- 
berg (both fc-space based) theories^. 

We consider the system described by the following set 
of parameters 62 — 61 = 2t,ti = t,t2 = 2t, ii2 = 0.05 
and n = 1.2, and start with the pairing interaction in 
the first band band only: Uu ^ 0, U12 = 0. Due to the 
weak hybridization ti2 — 0.05t there exist small coupling 
between bands. Figure ^ illustrates the changes in the 
order parameter Ai around intra-band (left panel) and 
inter-band (right panel) impurity. Note different pat- 
terns of changes in Ai. The intra-band impurity more 
strongly suppresses order parameter at the impurity site 
and leads to slight increase of it (with respect to the 
value for homogeneous system) at nearest neighbor sites. 
The intcr-band impurity scattering on the other hand di- 
minishes the order parameter at the impurity site and 
around it, but slightly less for next-nearest neighbors 
than for nearest-neighbors. In spite of its short range 
the inter-band impurity modifies the order parameter at 
distances larger than intra-band one. For the parameters 
used the clean system has Ai — 1.07i. At the impurity 
site one finds Ai(0,0) = 0.17< for V.mpin) - Klnp{r^) 
and Ai(0,0) = O.Ut for V^^pirl) = Vl,p{rl). 

In figure ([7]) we show the local quasiparticle density 
of states at the impurity site f = (0,0) and its nearest 
(0, 1) and next-nearest (1, 1) neighbor sites. The inter- 
band impurity induces states inside the gap. 

The interesting aspect of these studies is connected 
with the fact that the effect of V^^p depends on the sign 
of inter-band interaction Ui2. There is no similar de- 
pendence connected with intra-band impurities {Vj^p or 
V^^p). This is illustrated in the Fig. The changes 
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FIG. 6: The influence of intraband Vj^p{fi) (left panel) and 
interband Vil^p{fi) (right panel) impurity of the same strenght 
(=4t) on the local values of Ai(7^) in a superconductor with 
Uii = -3.5t, Ui2 = 0. 
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FIG. 7: The energy dependence of the local quasiparticle 
density of states projected onto orbital 1, at the impurity 
site (0,0) and its neighbours for system with intraband (left 
panel) and interband (right panel) impurities for Uii = —3.5t, 
Ul2 = 0. 



of the order parameters in the first and second bands 
clearly depend on the sign of inter-band interaction. For 
the clean system we have |A?| = 2.69t, |A§| = 2.07t. At 
the inter-band impurity site we have found |Ai| = 1.61t 
and IA2I = 1.21i for Uu = -5i, i.e. roughly 40% re- 
duction. On the other hand in the superconductor with 
U12 = +5t we find at the impurity site |Ai| = 0.20t 
and IA2I = 0.23t. The order parameters are suppressed 
few times stronger in the superconductor with repulsive 
inter-orbital interaction. 

Similar effect has earlier been noted^S within the weak 
coupling Eilenberger theory for finite (inter-orbital) im- 
purity concentration in the two band superconductors. 
Here we observe similar behavior already for the single 
impurity. The feature is farther discussed in the next 
section for a system with finite concentration of impuri- 
ties. Performing analytical studies of the Tc suppression 
in two band case the authorsiS have noted that strong 
suppression of superconductivity for finite concentration 
of inter-band impurities is to be expected for the inter- 
band couplings fulfilling the inequality 



Uu>- 



2iVi7V2 



(20) 



It other words much weaker suppression of is expected 
for attractive inter-band interaction (note that positive 
interactions are attractive in the notation of the paper 
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FIG. 8: The changes of the order parameters around inter- 
band impurity Vi]^p located at the center of the superconduct- 
ing cluster with Uu = — 3.5t. Left panels are for U12 = —5t, 
middle panels for Ui2 = 5t. The right panels show the maps 
of difference between values of the order parameter for attrac- 
tive and repulsive interband interaction. 




FIG. 9: The energy dependence of the quasiparticle density 
of states at the interband impurity Vimp{fi) = V^imp(?^) and 
the neighboring sites for attractive (left panel) and repulsive 
(right panel) interband pairing interaction. 



[39|). In more detail this is again illustrated in Fig. ([9]), 
which shows the energy dependence of the quasiparticle 
density of states at and near the impurity site for attrac- 
tive (left panel) and repulsive (right panel) inter-band 
interaction. In the later case the order parameter at the 
impurity site and around it is strongly suppressed and 
new states appear inside the gap. The intra-band impu- 
rity scattering (not shown) is not effective in suppressing 
superconductivity independently of the U12 sign. This is 
due to the same mechanism (Anderson theorem) as for 
single band s-wave superconductors^. 

It is also of interest to look at the relative phases of 
the order parameters near the impurity in the supercon- 
ductor with U12 > 0. As noted earlier for repulsive inter- 
orbital interactions the s±-wave pairing state^i^ is re- 
alised. This state has s-wave like order parameters on two 
Fermi surfaces, the phases of which differ by tt. It means 
that if the phase of Aa on one of the Fermi surfaces is 
(j), the phase on the other is (j> + tt. This phase relation 
is in fact responsible for strong suppression of supercon- 
ductivity by inter-band impurities. It turns out that at 
the V^l^p impurity site the phases of the gaps change by tt 
with respect to the phases in the bulk. Obviously, there 



7 



is no such effect for attractive inter-band interaction. 



V. MANY IMPURITIES IN THE TWO BAND 
SUPERCONDUCTOR 

In this section we consider the two band impure su- 
perconductor with inter-band and intra-band impurities. 
Our sample has a rectangular shape. It is L = 17 x 21 
sites large with a square lattice of unit lattice constant. 
It contains 20% inter-band or intra-band impurities ran- 
domly distributed. We are not averaging over distribu- 
tion of impurities, but rather calculate the local property 
at each site and present the result in the form of maps. 
To make the impurities more realistic we assume that 



they are extended V^^pifl) = V^^ fid with fid being a 
number from the Gaussian distribution at sites distance 
id ~ 1,a/2, 2 from the impurity. The superconductor 
studied in this Section is characterized by one active band 
with C/ii = — 3.5t and inter-band interaction \U\i\ = 5t. 
The other parameters are: 62 — ei = 2i, ti = t,t2 = 
2t,ti2 = 0.05 and n = 1.2. Energies are measured in 
units of t with respect to chemical potential. 

Figures pO)) and (jlip show the suppression of the or- 
der parameters in the two bands for intra-band and inter- 
band impurities, respectively. In all cases one observes 
similar patterns and large degree of (anti)correlation be- 
tween the impurity positions and the gap values. How- 
ever, the inter-band impurities suppress order parameters 
in both bands much less for U12 = — 5i than for opposite 
sign of this coupling. For negative U12 the phases of the 
order parameters on two Fermi surfaces are the same and 
the scattering of a pair from band 1 into band 2 is harm- 
less, as the superconductor as a whole looks like s-wave 
one, and is protected against impurities by the Anderson 
theorem'*'' . 

Even though the maps presenting suppression of the 
order parameters by the intra-band and inter-band im- 
purities shown in the Figures PH]) and pT|) look to large 
extend similar the big differences are observed in the lo- 
cal densities of states. They are shown in figures (fT2|) and 
([13]). Left panels of both figures show the local densities 
of states as function of energy along the line x = —7 in 
the first band. Middle panels present the results for the 
second band and the total DOS is plotted in right pan- 
els along the same cut. Intra-band impurities. Fig. 
induce large inhomogeneities, which show up as gaps in 
the local density of states of amplitude strongly changing 
from site to site. Sites close in space may have largely 
different gaps. Similarly, the inter-band impurities in a 
superconductor with large attractive inter-band scatter- 
ing also induce inhomogeneities. However, they are much 
smaller (recall for the same distribution and strength of 
impurities). The gaps also change from site to site, but 
more gradually i.e. on larger spatial scale. On the con- 
trary, the same inter-band impurities nearly completely 
destroy the superconductivity in the system with strong 
repulsive inter-band pair scattering i.e. in the s± state. 
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FIG. 10: Distribution of the intra-band impurity potential 
Vimpi'Ti) = Vi^pi^i) (left panel), local values of the order 
parameter in the first (middle panel) and second band (right 
panel) . The results do not depend if the inter-band scattering 
U12 = ±5t is of repulsive or attractive type. 
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FIG. 11: The maps of the local values of the order param- 
eters in the first band (left panels) and in the second band 
(right panels). The upper two figures are obtained for attrac- 
tive inter-orbital interaction U12 ~ —5t, while lower figures 
represent data for repulsive U12 ~ 5t interaction. 



Different reaction of the system with repulsive and at- 
tractive inter-band interactions to impurities may well 
be characterized by the dependence of the superconduct- 
ing transition temperature Tc or the average gap on the 
strength of Vimp. In the left panel of the Fig. we 
show the relative change of < Ai -I- A2 > in the im- 
pure system normalized to its clean value < Ai + A2 >o 
with increasing the strength of intra-band impurities 



on inter-band impurities Vi, 



imp 



^imp (dashed curve with triangles) and 
' ^imn ^ system with 



V\2 — — 5t (curve with full squares) or \J\2 = -l-5t (curve 
with dots). The right panel of that figure shows the 
changes in Tc with V^^^ for two signs of inter-band scat- 
tering. In a qualitative agreement with the results'^^ of 
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FIG. 12: Partial (left and middle panel) and total (right 
panel) density of states as a function of energy close to the 
Fermi energy at the sites along the line a; = — 7 of the same 
sample as shown in the figure (|10p . 
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FIG. 14: The left panel shows the dependence of the nor- 
malised average order parameter on the Vimp = 
(dashed curve with triangles) and on the Vimp ~ V^"^ with 
Uri = — 5t (curve with full squares) or U\2 ~ +5f (curve with 
dots). The right panel shows the changes in Tc with V^^ for 
two signs of inter-band scattering. 



Vi',™=2c.U„=-5l,x=-7 
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V|„,p=2t,U|2=-f5t, x=-7 
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FIG. 13: The energy dependence of the partial and total den- 
sities of states N(\,y,E) at the sites along the line x = —7 
of the sample shown in figure (|ll|l . The upper row corre- 
sponds to attractive inter-band interaction, while the lower 
to repulsive one. 



Kogan et al. we observe much stronger diminishing of 
Tc normahzed to its clean system value Tco for repulsive 
than for attractive Ui2. Similar dependence of Tc on im- 
purity strength for both signs of U12 at small disorder is 
attributed to our use of Bogolubov - de Gennes approach 
which is more suitable to study inhomogeneous systems 
than the Eilenberger or BCS theories. In BdG approach 
the condensate wave function may distort around the im- 
purity and this allows the system to keep its condensa- 
tion energy and Tc much higher than it would result from 
Abrikosov-Gorkov approach to impure superconductors, 
which does not allow for a local changes in the wave func- 
tion. 



VI. DISCUSSION AND CONCLUSIONS 

We have studied the model of two band superconduc- 
tor with intra-band and inter-band interactions. In the 
present paper we concentrated on some general aspects 
of the two band model. In particular we have found 
that the presence of Van Hove singularity in the density 
of states in one of the bands leads to strong enhance- 
ment of the superconducting transition temperature of 
intra-band only superconductor in the weak coupling, i.e. 
for Ao <C 1. However, neglecting changes of prefactor, 
the mere increase of effective interband coupling by Van 
Hove singularity seem to be not large enough to make 
the electron-phonon coupling responsible for supercon- 
ductivity, at least so in LaFeAsO for which the coupling 
constant matrix has been found^ with A12 — 0.093 and 
A21 = 0.124. Interestingly, our analysis suggests that for 
elevated values of Aq the effect of Van Hove singularity is 
to diminish Tc in comparison to the system without loga- 
rithmic enhancement of DOS. The model with both types 
of pairing interactions leading to s - or s±-wave symme- 
try displays number of features similar to that observed 
in real many band materials, particularly MgB2 and iron 
pnictides. 

MgB2 is a well established superconductor with two 
gaps^. It is believed to have one active band and rela- 
tively weak inter-band coupling. Its much slower than 
predicted by the Abrikosov-Gorkov theory^ suppres- 
sion of the superconducting transition temperature by 
impurities^i^ may point towards the inter-band char- 
acter of impurities and repulsive character of inter-band 
pair scattering {Ux2 > 0). Such a case is illustrated in 
the upper row of Fig. (fT5|) . 

Another prominent recent example of the many band 
superconductors is provided by the iron pnictides^. 
These superconductors seem to belong to different class 
of many band materials in which the inter-band interac- 
tion is dominant and the order parameter has different 
signs on different Fermi surface sheets^ - the s± state. 
The existing samples are certainly strongly disordered, 
as it can be inferred from large values of resistance just 
above Tc. Despite large disorder they are superconduct- 
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ing with quite large Tc- This suggests s- or s±-wave hke 
order parameter. As we have seen the superconductor 
with dominant attractive U12 interactions is quite robust 
against impurities, both of intra- and (and even more) 
inter-band type (c. /. figures (fT2)) and p3|) ). 

The robustness of the two band superconductors with 
attractive inter-band interactions to the impurities can 
be traced back to the Anderson theorem. On the other 
hand the s± state induced in two band model by repulsive 
inter-band interactions is characterized by large sensitiv- 
ity to inter-band impurity scattering. These results are 
in agreement with previous studies of similar models'^^!^. 

The maps plotted in the Figs. (|11I10[) are in qualitative 
agreement with recent scanning tunnelling microscopy 
(STM) studies for pnictide supcrconductors^^^i^. In 



these papers the relatively small variation of the local 
gaps have been observed with the average gap A = 
6 — ImeV and 2A(0)/fcsTc ~ 7 indicating strong cou- 
pling superconductivity. The detailed analysis of the 
STM spectra in pnictides will be the subject of future 
studies. 
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